A possible Efimov trimer state in a 3-component lithium-6 mixture 
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We consider the Efimov trimer theory as a possible framework to explain recently observed losses by inelastic 
three-body collisions in a three-hyperfine-component ultracold mixture of lithium 6. Within this framework, 
these losses would arise chiefly from the existence of an Efimov trimer bound state below the continuum of 
free triplets of atoms, and the loss maxima (at certain values of an applied magnetic field) would correspond to 
zero-energy resonances where the trimer dissociates into three free atoms. Our results show that such a trimer 
state is indeed possible given the two-body scattering lengths in the three-component lithium mixture, and gives 
rise to two zero-energy resonances. The locations of these resonances appear to be consistent with observed 
losses. 



The experimental realisation of ultracold Fermi gases 
have let us explore many fundamental aspects of few and 
many-body physics. In particular, the study of mixtures 
of fermionic atoms in two different spin components have 
led to the observation of superfluid paired phases such 
as molecular Bose-Einstein condensates (BEC), Bardeen- 
Cooper-Schrieffer (BCS) superfluids, and their crossover 
im 121 [3] [4] [51 . Recently there has been some theoretical in- 
terest in Fermi systems with three different spin components 
E H 121 Uni [TT|, which can present analogies with colour 
superfluidity in QCD [8 |. Recent experiments 1 121 [13] have 
been performed with ultracold mixtures of lithium 6 atoms 
prepared in the lowest three hyperfine states |2), and 
|3). They indicated that when an external magnetic field is 
applied, strong losses due to three-body inelastic collisions 
occur over a wide range of magnetic field intensities. On 
the other hand, such losses are not observed when only two 
of the three hyperfine components are mixed. Therefore the 
observed inelastic collisions are related to the specific scat- 
tering channel involving three atoms in the three different 
hyperfine components. The magnitude of the inelastic col- 
lisions is characterised by a loss rate coefficient ,J^, defined 
by the rate equation 

= -XniYijU],, for {i,i,k) = (1,2,3) 

where n\, n^, and ht, are the densities of each kind of atoms. 
The variation of the measured loss rate coefficient with re- 
spect to the intensity B of the applied magnetic field is shown 
in the bottom panel of Fig. [T] It reveals a peak around 
B =130 G which suggests an enhancement due to a reso- 
nance of three colliding atoms with a three-body bound state. 
We can envisage two kinds of resonance, depending on the 
origin of such a three-body bound state. 

Firstly, the three-body bound state may originate from an- 
other hyperfine channel and couple by hyperfine interaction 
to the three-body scattering continuum in the hyperfine chan- 
nel |1)|2)|3). Since the bound state and scattering state be- 
long to different hyperfine channels, they have different mag- 
netic moments and become resonant only around a particular 
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intensity of the magnetic field which brings them to the same 
energy. This situation would correspond to a "three-body 
Feshbach resonance" |14|, a generalisation of the now well- 
known two-body Feshbach resonances which occur for two 
scattering atoms around certain magnetic field values ifTSll . 

As a matter of fact, wide two-body Feshbach resonances 
are present in lithium 6 over the range of magnetic field in- 
tensities where the 3-body losses are observed. As a result, 
the two-body scattering lengths between two atoms in differ- 
ent states, namely |1)|2), |1)|3), and |2)|3), are modified by 
the applied magnetic field. The dependence of these three 
scattering lengths on the magnetic field intensity is shown 
in the top panel of Fig. [T] Because of this dependence, the 
interactions between three atoms colliding in the hyperfine 
channel 1 1) |2) |3) are also modified by the magnetic field, and 
it may happen that a 3-body bound state supported by these 
interactions within the same hyperfine channel is brought to 
the threshold of its three-body scattering continuum at a cer- 
tain magnetic field value, causing a "shape resonance". This 
constitutes the second possible kind of resonance. 

Studying both kinds of resonance theoretically is involved 
and requires an extremely accurate knowledge of the inter- 
actions between atoms. However, in the case where the two- 
body scattering lengths are much larger than the range of the 
interatomic interactions, it is possible to predict the struc- 
ture of the three-body bound states near threshold and their 
"shape resonance" simply in terms of the scattering lengths 
and a short-range 3-body parameter. This was pointed out by 
V. N. Efimov fTSl, and the corresponding three-body bound 
states, known as "Efimov trimers", are purely quantum- 
mechanical states which enjoy special properties such as dis- 
crete scale invaiiance as the scattering lengths are varied. In 
particular, their energy spectrum forms an infinite series with 
a point of accumulation just below the continuum threshold 
when the scattering lengths become infinite. So far, signa- 
tures of Efimov trimers of identical bosons have been ob- 
served in helium 4 1 19|, cesium 133 li20]|2TI . and potassium 
39 [22] and have been assigned in the first two cases to the 
ground state of the Efimov series lfT8ll24l . and to the ground 
and first excited state in the case of potassium. Similar signa- 
tures of heteronuclear bosonic Efimov trimers of potassium 
and rubidium were recently reported in Ref. 1231 . while the 
evidence of fermionic Efimov trimers is yet to be found. It 
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Figure 1 : Top panel: variation of the two-body scattering lengths 
ai2, fli3, and a23 for the lowest three hyperfine components of 
hthium 6 as a function of magnetic field. These curves were cal- 
culated by P. S. Julienne and taken from Ref 1 12|. Middle panel: 
Energy of the Efimov trimer (solid curve) just below the three-body 
threshold as a function of magnetic field with Ay = (0.42 flo) 
The shaded area corresponds to the width of the Efimov state, i.e. 
the imaginary part of its energy, for rj = 0.115. The estimated en- 
ergy of possible resonant trimers from all other spin channels with 
the same total projection m/r = — 3/2 is indicated by dashed curves. 
Bottom panel: Experimental (dots, taken from Ref. [12|) and the- 
oretical (curves) three-body inelastic collision loss rate coefficient 
as a function of magnetic field. The dashed curve is obtained by 
adjusting the short-range loss parameter rj to fit the experimental 
data (rj = 0. 157), and the solid curve by adjusting it to fit the shape 
of the left peak (t] = 0. 1 15). 



was suggested by the authors of Refs. lfT2l[T3l that their ob- 
servations in three-component lithium 6 might be the man- 
ifestation of an Efimov trimer of distinguishable fermions. 
The purpose of this Letter is to test this hypothesis using the 
Efimov theory with the parameters of the experiments. 

First, it should be noted that the two-body scattering 
lengths in the experimental conditions are indeed quite large, 
from about -100 to -1000 oq (where oq 5.292 10"" m is 
the Bohr radius), but not always much larger than the range 
of the atomic interactions, typically given by the van der 

Waals length €,,c/w = (wQ/^)^''^ ~ 60 qq, where Q is the 
van der Waals dispersive coefficient, m is the mass of lithium 
6, and h is the reduced Planck's constant. Therefore, the 
applicability of Efimov theory is questionable, especially at 
low (around 100 G) and high (around 500 G) magnetic field 
values where one or two of the scattering lengths become 
small. However, in the intermediate region, the necessary 
conditions for the existence of Efimov trimers are met. 

The details of the Efimov theory can be found in Refs. lfT6l 
[TSl . It essentially treats the free three-body problem with 



boundary conditions at short distance imposing the known 
two-body scattering lengths between each pair of atoms. The 
three-body wavefunction ^'{R,a, 9) is expressed in terms of 
the hyperradius R of the three-body system (a measure of 
the global distance between the three atoms) and the two 
angles a and 9 describing the geometrical configuration of 
the atoms. More precisely, if we denote by r the relative 
vector between atom 1 and 2, and by p the relative vector 
between atom 3 and the centre of mass of atoms 1 and 2, 
then = + |p^, 9 is the angle between r and p, and 
tana = \/3r/(2p). Using the Faddeev decomposition ifTTI . 
and restricting ourselves to the case of zero total angular mo- 
mentum which is the most favourable for the Efimov effect 
to occur 1161, the wave function can be written as 



R^ V sin2a+ 



sin2a_ 



sin 2a 
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where sin2a± = [l — (jcosai ^sinacos0)^] ' 

Note that the functions j''' depend on only one hyperan- 
gle, because we assumed that the total angular momentum is 
zero. They satisfy the free Schrodinger equations 
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for a given energy E, with the boundary conditions 
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where /, j, k is any permutation of 1 , 2, 3. 

The first boundary condition Q imposes the form of the 
wave function in the two-body sectors consistent with the 
known two-body scattering lengths an, flis, and 023 be- 
tween the three kinds of atoms. The last boundary condi- 
tion Q fixes the logarithmic derivative of the wave func- 
tion at short hyperradius Rq ^ Oij to some value A(7?o) in- 
dependent of the energy E. In this region, the last term 
of Eq. Q is negligible, and one can show that the so- 
lution of Eqs. ilt takes the separable form x^^\R,a) « 
sin(|so|ln(/r7J)+A)sin(so(7r/2-a)), where \fmElh, 
A is a phase shift, and so ~ 1.00624/ is the imaginary solu- 
tion of the equation — socos(so7r/2) +8/\/3sin(so7f/6) = 0, 
which follows from applying the boundary condition Eq. (|2]l. 
In order for A to be energy independent, the phase shift A 
must be of the form — |so| ln(-^^/Ao) where Aq is some fixed 
wave number. This sets 



A(i?o) 



Nl 



cot[|so|log(i?oAo)]. 
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This way, the choice of Rq is arbitrary. The only requirement 
is that Rq should be much smaller than the scattering lengths. 
In our calculation, we fixed Rq to 1 qq. Thus, the only free 
parameter of the theory is Aq. It captures the effects of the 
unknown short-range three-body physics on the wave func- 
tion at larger hyperradii. 

We solve these equations by discretising the arguments 
{R,cc) of the functions x'"'^ on a 2-dimensional grid, eval- 
uating the derivatives by finite differences, and diagonalising 
the resulting matrix corresponding to the Hamiltonian of the 
system. The maximal value of R is set to 50,000 oq, which 
is on the order of the particle spacing in the experiments. At 
each value of the magnetic field, the three scattering lengths 
ai2, ai3, and a23 are obtained from the top panel of Fig. [T] 
The only unknown quantity is Aq. First, we make the as- 
sumption that it does not depend on the magnetic field. In 
reality, it may actually depend on it, but provided that no 
accidental resonance occurs, it is reasonable to assume that 
its variations are less pronounced than that of the scattering 
lengths, at least in first approximation. Second, we make the 
fundamental assumption that the observed three-body losses 
are due to a shape resonance with an Efimov trimer The 
measured three-body loss rate coefficient as a function of 
the intensity of the magnetic field shows a distinctive peak 
around B = 130 G - see the bottom panel of Fig.[T] Assuming 
that this is the point where an Efimov state reaches the con- 
tinuum threshold (its binding energy goes to zero), we find 
that we should adjust Aq to about (0.42 ao) Once Ao is 
fixed, we can obtain the eigenstates and their energy around 
the threshold, namely the Efimov trimers and the three-body 
scattering states. 

In the middle panel of Fig. [TJwe plotted the energy of the 
Efimov trimer just below the continuum threshold as a func- 
tion of magnetic field. By construction, the trimer appears at 
B = 130 G. Interestingly, its binding energy increases until 
B — 350 G, and then decreases until the trimer reaches the 
continuum again, causing a second zero-energy resonance 
around B — 500 G. This simply results from the magnetic 
field dependence of the scattering lengths. 

Since the two-body scattering lengths are all negative in 
this range of magnetic field, there is no two-body bound 
states just below the two-body continuum. As a result, when 
inelastic three-body collisions occur, two of the three atoms 
have to form a deeply bound dimer A direct calculation of 
the rate for such processes would require a detailed analysis 
of the deeply bound dimers. However, we can easily calcu- 
late its enhancement by the Efimov resonance. Indeed, by 
assigning a complex value |Ao|e"' to the 3-body parameter 
Ao, and therefore a complex value to the logarithmic deriva- 
tive A, one can impose a probability loss at short hyperra- 
dius (R = Rq) in order to model the overall effect of losses 
by recombination to deeply bound dimers |18|. This makes 
any scattering state ^ quasistationary, and by calculating the 
time variation of the total probabiUty / \^\^d^Rd^r, one can 
derive the following 3-body loss rate coefficient: 

2h 

JT = — |ImA| X p{Ro). (5) 
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Figure 2: Probability density p{R) - defined by Eq. |6]( - of the low- 
est three-body scattering state for different values of the magnetic 
field. In all cases, the wave function is normalised to be asymp- 
totically equal to the noninteracting limit 8J2{KR)/ (KR)^, where 
J2 is the Bessel function, and K is the wave number. This limit is 
indicated by the dashed curve. Here, the wave number K is set by 
the size of the numerical grid. 

As one would expect, this coefficient is proportional to the 
imposed velocity Y — ^|ImA| at R = Rq and the probabil- 
ity density p{Rq) of finding three atoms at that hyperradius, 
defined by the hyperangular average 

p{R) = sh^R^J^'^l^ sinlafda J\m9de\^'{R,a,9)\^, 

(6) 

where ^ is normalised to be asymptotically equal to the non- 
interacting limit %J2{KR) / {KR)~ , J2 being the Bessel func- 
tion. This probability density, and thus the inelastic pro- 
cesses, is strongly increased at short distance by the presence 
of an Efimov trimer just below threshold. Physically, this is 
due to the fact that the three atoms almost bind during their 
collision, and therefore spend more time together. One can 
see in Fig. |2] that while p{R) is unaffected at large distance, 
it changes significantly at short distance when the magnetic 
field is varied around the zero-energy resonance. 

The calculated loss rate coefficient of Eq. Q is plotted in 
the bottom panel of Fig. [T] We obtain a profile delimited by 
two peaks. In between the two peaks, we observe a plateau 
due to the presence of the Efimov state below threshold. Out- 
side, the probability becomes very low, due to the absence of 
a near-threshold Efimov trimer. From the energy spectrum, 
it is clear that the second peak around 500 G is related to 
the second zero-energy resonance with the Efimov state. We 
can then adjust the value of Tj to fit the experimental data. 
It should be noted that Tj has two effects: it sets the overall 
magnitude of the loss rate coefficient, as can be seen from 
(|5]l, and it also smoothes the peaks (because they are natu- 
rally broadened by the losses). It turns out that there exists 
a range of values for T] which can approximately fit both the 
overall magnitude of the rate coefficient and its shape. Us- 
ing a least-square minimization method, we found that we 
should set tj to 0. 157 in order to obtain the best fit to the ex- 
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perimental data. Fitting only the shape of the first peak with 
that of the experimental peak at B = 130 G, we obtain the 
best agreement for rj =0.115. 

The behaviour of the calculated three-body decay rate co- 
efficient is very reminiscent of the measured one, which has 
a similar profile between 130 and 500 G. This suggests that 
the local maximum near B = 500 G found in the experiment 
is caused by a second resonance. The agreement with the 
observations is only approximate however, as the experimen- 
tal data show a much more diffuse local maximum. As we 
noted earlier, the Efimov theory is not strictly applicable to 
the present system, and it is expected that short-range cor- 
rections are needed for a better agreement. A magnetic-field 
dependence of Aq might also play a role. These effects may 
very well explain the remaining discrepancies. Yet, it is quite 
remarkable that the Efimov theory already provides a basic 
description which qualitatively explains the experimental ob- 
servations. In this interpretation, a relatively "pure" Efimov 
trimer state is expected around 300 G (where the scatter- 
ing lengths are largest) and connects continuously to trimer 
states which are partly affected by short-range physics ("im- 
pure" Efimov trimer) and eventually dissociate in the contin- 
uum at both resonances. We find that the pure Efimov trimer 
at B =300 G has a binding energy of about 2nh x 10 MHz and 
a width (imaginary part of the energy) of about 2nhx3 MHz, 
corresponding to a lifetime of about 50 ns. However, these 
values result from an adjustment near B = 130 G, where the 
Efimov theory may need corrections. Thus, the actual values 
might be slightly different. 

The present interpretation may be confirmed or refuted by 
further experimental investigation, in particular direct obser- 
vation of the trimers below threshold - for example, by ra- 



diofrequency spectroscopy. For comparison with the first 
interpretation discussed in the introduction (Feshbach res- 
onance with another hyperfine channel), we calculated the 
expected energy of possible resonant trimers coming from 
other channels, assuming that their magnetic moment is sim- 
ply the sum of the magnetic moments of the three separated 
atoms. Since the interaction conserves the projection of the 
total spin, we considered only the channels with the same 
projection mp = —3/2 as the incoming channel |1)|2)|3), 
and shifted the energy of the possible resonant trimers such 
that it crosses the 3-body threshold of the incoming chan- 
nel at B = 130 G. The resulting energies as a function of 
magnetic field are plotted in Fig. [T] as dashed curves. One 
can see that they have a monotonic behaviour, and depart 
steeply from the continnum threshold. On the other hand, if 
the Efimov interpretation is correct, the trimer is expected to 
follow a different and rather unusual behaviour: it connects 
to the continuum via two zero-energy resonances, resulting 
in an energy minimum as function of magnetic field. Thus 
we have provided a significant difference between the two 
scenarios, which may serve as a test in future experiments. 

We are grateful to S. Jochim for providing the experimen- 
tal data shown in Fig. 1 . 

Note added: Upon finishing this Letter, we became aware 
of the work by E. Braaten et al. Il25l . which uses the 
Skorniakov-Ter-Martirosian equations to directly calculate 
the loss rate from scattering amplitudes. After this Letter was 
submitted, similar results were reported in Ref. [26 1, using a 
different model involving the formation of a trimer Both of 
these works yield results which are consistent with the loss 
rate coefficient shown in Fig. 1. 
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